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Abstract 

We show that a model of the hippocampus introduced 
recently by Scarpetta, Zhaoping & Hertz ([2002] Neural 
Computation 14(10):2371-96), explains the theta phase 
precession phenomena. In our model, the theta phase 
precession comes out as a consequence of the associative- 
memory-like network dynamics, i.e. the network's ability 
to imprint and recall oscillatory patterns, coded both by 
phases and amplitudes of oscillation. The learning rule 
used to imprint the oscillatory states is a natural gener- 
alization of that used for static patterns in the Hopfield 
model, and is based on the spike time dependent synaptic 
plasticity (STDP), experimentally observed. 

In agreement with experimental findings, the place cell's 
activity appears at consistently earlier phases of subse- 
quent cycles of the ongoing theta rhythm during a pass 
through the place field, while the oscillation amplitude 
of the place cell's firing rate increases as the animal ap- 
proaches the center of the place field and decreases as 
the animal leaves the center. The total phase precession 
of the place cell is lower than 360°, in agreement with 
experiments. As the animal enters a receptive field the 
place cell's activity comes slightly less than 180° after the 
phase of maximal pyramidal cell population activity, in 
agreement with the findings of Skaggs et al (1996). Our 
model predicts that the theta phase is much better corre- 
lated with location than with time spent in the receptive 
field. Finally, in agreement with the recent experimen- 
tal findings of Zugaro et al (2005), our model predicts 
that theta phase precession persists after transient intra- 
hippocampal perturbation. 



1 Introduction 

In this paper we will show how the model of the hippocam- 
pus introduced in Scarpetta et al. (2002) is able to ex- 
plain the theta phase precession phenomena: the observa- 
tion that oscillation phases of the activities of hippocam- 
pal place cells systematically shift relative to the theta 
rhythm as the animal traverses the place fields (O'Keefe 
Recce 1993; Skaggs et al 1996). 

The first part of the paper reviews the model, its learn- 



ing rule, and the recall dynamics. Briefly, the model's 
learned states are encoded by both the amplitudes and 
the phases of the oscillating neural populations, enabling 
more efficient and robust information coding than in con- 
ventional models of associative memory. The learning rule 
to imprint such oscillatory states is a natural generaliza- 
tion of that used for static patterns in the Hopfield model, 
and is based on the STDP synaptic plasticity observed 
experimentally. In particular, long-term potentiation and 
depression of the synaptic efficacies depends on the rela- 
tive timing of the pre- and post- synaptic activities during 
learning. A further step is introduced to allow imprinting 
of correlated patterns properly. 

The second part of this paper shows how the model 
accounts for theta phase precession observed in rat hip- 
pocampus. During learning of each spatial location, our 
model imprints an oscillatory pattern of specific ampli- 
tude and phase relationship between neurons by modify- 
ing neural synaptic connection strengths. After learning, 
the network's response to external inputs reveals the en- 
coded phase information, and thus phase precession, even 
though the external inputs do not have the phase infor- 
mation. The phase precession effect comes out from the 
crucial ability of our model to store both the amplitudes 
relationship and the phases relationship of the oscillatory 
patterns in the synaptic connections. Both the bell-shaped 
firing tuning curves (receptive fields) and the phase pre- 
cession phenomena are a result of the network synaptic 
connections. In particular, the phase precession is an ef- 
fect of the overlap between the imprinted place fields: if 
adjacent receptive fields do not overlap the firing tuning 
curve of each place cell is preserved but the phase preces- 
sion does not occur. 

The model agrees with the observation (Skaggs et al 
1996) that the first spikes, as the animal enters the recep- 
tive field, come near the phase of the minimal pyramidal 
population average activity. The total phase precession of 
the place cells can be close to 360°, but not larger than 
360° as found by O'Keefe Recce (1993) and Skaggs et al 
(1996). In our model, the theta phase of the place cell 
is a function of the animal location, almost insensitive to 
the animal' running speed (for speeds that do not exceed 
a critical value related to the network connections and pa- 
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rameters), and therefore the theta phase shows a better 
correlation with the animal location than with time, in 
agreement with the experimental results of O'Keefe Recce 
(1993) and Huxter et al. (2003). 

Numerous computational models (Burgess et al 1994, 
Tsodyks et al 1996, Wallenstein et al 1997, Jensen and 
Lisman 1996, Bose et al 2000, Booth et al 2001, Yam- 
aguchi 2003, and Lengyel et al 2003 ) have attempted to 
account for the mechanisms that underlie the phase pre- 
cession observation, but a fully-satisfactory explanation is 
still lacking. Recently, to contrast the predictions of differ- 
ent models of phase precession, Zugaro et al (2005) have 
transiently turned off neuronal discharges for up to 250 ms 
and reset the phase of theta oscillations by stimulating the 
commissural pathway in rats. After recovery from silence, 
phase precession continued: the theta phase of the place 
cell immediately after the recovery was correlated with the 
new location of the animal, despite the transient silence 
and the theta phase reset. Our model is in agreement with 
these recent experimental findings. 

2 The model 

The model structure and elements are based on the phys- 
iological and anatomical findings in the CAS hippocam- 
pal region. These regions contain the principle pyramidal 
cells and the inhibitory interneurons. The pyramidal cells 
project long range axons to other pyramidal cells and in- 
terneurons, whereas the interneurons project more locally. 
The model is based on the one introduced and analyzed 
by one of us with Li Zhaoping and John Hertz in a recent 
paper (Scarpetta et al 2002). We review here its main 
features, more details can be foimd in the paper above 
mentioned. The state variables, modeling the membrane 
potentials, are u = {wi,..., w„} and v = {v\, . . . ,Vn} 
respectively for the excitatory and inhibitory imits. (We 
denote vectors by bold fonts.) The unit outputs, repre- 
senting the probabilities of the cells firing (or instanta- 
neous firing rates) are given by (/„(ui),..., gu{un) and 
5^(111), ... , gv{vn)i where g„ andg^ are sigmoidal activa- 
tion functions that model the neuronal input-output rela- 
tions. The equations of motion are 

iii = -aui - P^Qvivi) + ^ Jijguiuj) + h, 

Vi = -ai^i +7»V(wi) +^ W°.Sfu(wj). (1) 

where is the membrane time constant (for simplicity 
assumed the same for excitatory and inhibitory units), J^j 
is the synaptic strength from excitatory unit j to excita- 
tory unit i, W^j is the synaptic strength from excitatory 
unit j to inhibitory unit i, jS^ and 7? are the local in- 
hibitory and excitatory connections within the E-I pair i, 
and Ii{t) is the net input from other parts of the brain. 
The model units represent local populations of biological 
neurons that share common input, so the number of neu- 
rons represented by an excitatory unit may be different 
from the number of neurons represented by an inhibitory 



unit. In this minimal model we omit inhibitory connec- 
tions between pairs, since the real anatomical long-range 
connections appear to come predominantly from excita- 
tory cells. The sensory input Ii{t) to the system, which 
drives the excitatory units, has a static part li and an 
oscillatory part SIi{t) modulated at the theta frequency. 
The static part I of the input determines a fixed point 
(u,v), given by the solution of equations u = 0, v = 
with 1 = 1. 

Linearizing the equations (1) around the fixed point 
leads to 

iii = -aUi - piVi + ^ JijUj + Sli 

j 

Vi = -avi + jiUi + ^ WijUj (2) 
i 

where Ui and Vi are now deviations from the fixed point, 
(51 = I — I is the oscillatory part of the sensory input, 
A = .g;(«0/9°, 7. = g'uiu^)l°, W,, = g'^{u,)W^^, Ji^ = 
g'u{uj) Jij. Henceforth, for simplicity, we assume /3i = /3, 
7j = 7, independent of i. 

Eliminating the Vi from (2), we have the second order 
differential equations 

ii -I- (2a - J)u + [q2 -aJ+ /3(7 + W)]u = {dt + a)SI. (3) 

(We use sans serif to denote matrices) or, equivalently, 

[{dt + af + /37] u = Mu + {dt + a)51 (4) 

where M = (5* a) J - /3W , (5) 

The terms in the square bracket describe the local E- 
I pair dynamics, while M gives an effective coupling be- 
tween the oscillating E-I pairs. Learning imprints pat- 
terns into M through the long range connections J and W. 
After learning, u depends on how I decomposes into the 
eigenvectors of M. Thus the network can selectively am- 
plify or distort I in an imprinted-pattern-specific manner 
and thereby functions as input mapping into the subspace 
spanned by the imprinted patterns. 

Following Scarpetta et al (2002), we distinguish a learn- 
ing mode, in which the oscillatory patterns are imprinted 
while the connections J and W are ineffective, from a re- 
call mode, in which connection strengths are effective and 
fixed. With appropriate learning kernel, the network dur- 
ing recall responds strongly (resonantly) to inputs simi- 
lar to those learned or linear combinations of them, and 
weakly to other inputs not related to the imprinted pat- 
terns. The phases and the amplitudes of the oscillatory 
response will be a linear combination of the phases and 
the amplitudes, respectively, of the imprinted oscillations 
recalled by the input. 

2.1 Learning oscillatory patterns 

First, let's imprint a single oscillatory pattern, P'^{t) = 
^A^g-i^o* _|_ where c.c. denote complex conjugate, and 
is a complex vector with components = |^j'|e*'^J . 
Following Scarpetta et al (2002), the synaptic strength 
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Jij is shaped by the correlations of the activities, Ujit') 
and Ui{t), in pre-synaptic and post-synaptic cells through 

-^'^ = (6) 

where Aj{t — t') is the learning kernel that measures the 
strength of the synaptic change at time delays t = t — t' 
between the pre-synaptic and post-synaptic activities. In 
the case when Aj{t — t') = S{t — t') the learning rule 
becomes the conventional Hcbbian learning with oc 
J dtui{t)uj(t) used in Li and Hertz (2000). To model the 
experimental results of STDP (Markram et al 1997, Bi & 
Poo 1998) the kernel A{t) was taken as an asymmetric 
function of r, mainly positive (LTP) for r > and mainly 
negative (LTD) for r < 0. Analogously for the synaptic 
strength Wij we write 
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dt 



dt' Vi{t)Aw{t - t')uj{t'). (7) 



Since the connections J and W arc ineffective during learn- 
ing, the responses Ui and Vi are proportional to the im- 
printing input Pi{t). Substituting them into eqs. (6) and 
(7) yields connections 



-Re[Aj{u;o)inr 



27^ 



a — iwo 



(8) 



where 

Aj^w{i^) 
with xo(w) 



/oo 
dTAj,w{T)e-'^^ (9) 
-oo 

a — iwo 



(10) 



Note that ImAj^vy(w) = if Aj^w{t) is symmetric in r 
and that He Aj^w{'^) = if Aj,w{t) is antisymmetric. 
The dependence of the neural connections J and W and 
the oscillator couplings M on is just a natural gener- 

alization of the Hebb-Hopficld factor for (real) static 
patterns. This becomes particularly clear for kernels sat- 
isfying the following matching condition: 



Aj{oJo) 



+ oj: 



;Aw{ijJo), at IV = LOq. (11) 



Then the oscillator coupling simplifies into a familiar 
outer-product form for complex vectors ^: 



-2h;oAj{u;o)^^^f*/N, 



(12) 



In the following we will always consider the case of 
matching kernels (11). 

We are interested in generalizing eq. (8) and (12) to 
the case in which multiple correlated patterns have to be 
imprinted. To construct the corresponding matrices for 
multiple patterns, /x = 1,..,N, we use the rule pro- 
posed by Personnaz et al (1985) and Diederich and Opper 



(1997) (see also Hertz el al 1991) for learning static corre- 
lated patterns without errors in the Hopficld model. Thus 
we write, for N, linearly independent, correlated patterns 



J,,^Re (^i(u;o)C^(g-^)M.?r) 
where the N x N matrix Q, defined by 



(13) 



(14) 



is the natural generalization of the correlation matrix Q 

used by Personnaz et al (1985) and Diederich and Opper 
(1997) for the Hopfield model. The rule (13) is not local 
and not iterative. However it has been proved in Diederich 
and Opper (1997) that, in the Hopfield framework, there 
exists a local rule that uses only successive presentation 
of one pattern at a time that converges exactly to the 
synaptic efficiencies obtained with the Q-rule. Work is 
in progress to generalize their procedure, and obtain the 
same result also in our case. In the following, we have used 
the rule (13). With the rule (13), we imprint N oscillatory 
patterns = ^''e"^"*, all with the same frequency coq. 
Then the M matrix becomes: 



Mi. 



(15) 



2.2 The recall mode 



In the recall mode the network dynamics is governed by 
the matrix M, learned during the learning mode. Keep- 
ing in mind cq. (15), it is clear that the imprinted oscilla- 
tions (even though not-orthogonal) are eigenvectors of the 
connection matrix M, 



-2icjoi(t^o)r 



(16) 



Thus, after the transient (which is governed by the net- 
work reaction time that depends both on the membrane 
time constant of the neurons, a, and the network synap- 
tic connections) , the response u to an input that perfectly 
match one of the imprinted vector is 

u = xl (17) 
with the linear response coeSicient (or susceptibility) 



X = 



a — ilvq 



a"^ + Pj — u)q — 2i(Joa + 2iu)oA{uo) 



(18) 



A properly chosen learning kernel A{t) satisfying 

a'^+f3^-LJ^ + 2ia;oRe [A{ljo)] 
a - Im [^(wo)] — >■ 0+ 



(19) 



would produce a very strong susceptibility x and thus a 
resonant response to match input. Meanwhile, an input 
pattern I outside the subspace of the imprinted vectors 
evoke a weaker response u = xd- In general, an input 
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/ = ^e"^* that overlaps with several imprinted patterns 
will evoke a correspondingly mixed resonant response 

u = Xl|| + XoI± (20) 

where I = I|| + I±, I|| = E^(»7^IO?^e-''"* = 

N-\Ej vr^j)^"^-'^'' and j,^ = Q.^e- 
This feature enables the system to interpolate between 
imprinted patterns, and to perform an elementary form of 
generalization from the learned set of patterns. As it will 
be shown in the next section, the phase precession as the 
animal moves from one location to the other is related to 
this interpolation mechanism. 

An example of learning kernel A{t) that satisfy the res- 
onance conditions (19) at frequency t^o = 10 Hz (theta 
range) is shown in Fig. 1. This learning kernel is able 
to imprint oscillations with frequency luq = 10 Hz, and to 
make the network resonant at theta frequencies. It could 
be of interest to note that the kernel A{t) shown in Fig. 
1 satisfies the resonance condition also at ojq = 40Hz, i.e. 
in the gamma range. However we will not exploit this last 
property in this paper. 

3 Imprinting locations and the 
theta phase precession 

Let's consider N locations, arranged in a 1-dimensional 
space with periodic boundary conditions (such as in a cir- 
cular or triangular track). For each of these N locations 
x^, ^ = 1,..,A^ we imprint, in the synaptic strengths of 
the network, an oscillating pattern P^, with specific ampli- 
tude and phase relationship between its components. Us- 
ing eq. (13) in order to imprint N (linearly independent) 
patterns, one need a number of network units n > N. In 
the following we use a network with n=2N excitatory units 
and the same number of inhibitory units, that follows the 
equation of motion (1). For each location an oscillat- 
ing pattern P'' , whose amplitude profile is peaked in the 
excitatory unit 2/x, is imprinted in the network. In par- 
ticular, the imprinted pattern associated with location x^ 
has components Pf (t) = ^^^e"''^"* + c.c, fi = 1,2,. .,7V, 
j = 1,2,..,2N, where = A(2/^ - j)e*'^(2p-j) . the am- 
plitudes A and the phases </> of the pattern ^^fj, on each 
unit j are functions of the distance between the unit j 
and the unit 2/i that is associated with location .x^. In 
particular, A() is an even function decreasing quickly with 
distance, while ^{) is an odd function. This means, when 
imprinting a location x^, only units located close to unit 
2/x are excited by the imprinting inputs. In such a man- 
ner, after the learning process, the N excitatory units U2/i, 
that we will call place cells, will be associated with a re- 
ceptive field centered at x^, for /x = 1, 2, . . . , A^. The other 
N units (with odd index, called auxiliary units) are neces- 
sary in order to allow the network to store the N desired 
patterns properly. In all simulations we used A(0) = 1, 
A(±l) = 0.3 A(±2) = 0.3, A(x) = otherwise, and 
0(0) = 0, (t){±l) = T'^Arad and (/)(±2) = =F2.5md. The 
network has a quite robust behavior with respect to the 
choice of the functions A and ^. 



In our model the external inputs in the recall mode 

resembles but are not exactly the imprinted patterns. 
When the animal is exactly at location x^ , the input pat- 
tern 51 is given by = L'^e""^"* where L'' is a vector 
with components = 5j^2fi- This pattern conveys the 
animal location in the amplitude modulation but has no 
phase coded information. When the animal is at position 
X between location and location x^+i, x^ < x < x^+i, 
we take the input vector as linear combination of L'^ and 

L- = (l-|a;-x^|)L'' + |x-x^|L''+i (21) 

L is then normalized. Thus, if the animal moves from 
location x^ at time i = to the adjacent one x^+i with 
speed v{t), the current location will be a function of time 
x{t) = Xf^i + Jj* v{t)dt, and the input vector will be ex- 
pressed by eq. (21) where x is replaced by x{t). The 
vector L is always normalized. Fig. 2B shows the input 
amplitude Lj (t) when animal moves along the track with 
constant speed. 

A animal at location x receives an input pattern 
that resembles imprinted pattern for x^ w x. Hence, 
the network response to input would amplify the im- 
printed pattern that resembles most (see previous 
section) , and thus exhibits in the responses the non-trivial 
phase relationship among units that was stored in the cou- 
plings. When animal is in between two stored locations, 
the response will be a combination of the two imprinted 
patterns (see eq. (20)), and in particular, the phase of 
the response will be a combination of the phases of the 
two imprinted overlapping patterns which depends from 
the animal position. So, even though the input carries in- 
formation about current location only in the amplitudes 
modulation, the associative memory nature of our model, 
having phase-coded locations imprinted in the synaptic 
couplings, allows the network to show both amplitude and 
phase coded informations in the response. 

Since the theta rhythm is the local field potential, in 
our model it is simply the mean field. We compute the 
theta rhythm doing the average of all the excitatory units 
activities M.F. = 1/2 A ^^^^ Uj{t). This choice allows us 
to compare the the absolute phase of theta at which place 
cells fire during the place field traversal, with the experi- 
mental findings of Skaggs et al. (1996). Indeed, the phase 
zero of the theta rhythm in the experimental work of Sk- 
aggs et al 1996 is defined to be the point in the theta cycle 
corresponding to the maximal pyramidal cells population 
activity (i.e. the average over the entire dataset of pyra- 
midal cells), and the EEG is opportunely shifted in order 
to agree with this definition of theta rhythm phase. 

4 Simulation Results and compar- 
isons with experimental findings 

In all simulations reported here we used the kernel A(t) 
shown in fig. 1 that satisfy the constraints (19) when 
loq = lOHz, i.e. when the frequency is in the theta range. 

After having imprinted A = 10 locations in the net 
with n = 20 excitatory units, we simulate the network 
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dynamics while the animal moves continuously in space. 
Fig. 2 shows in red the activity of the two place cell units 
{u2ij, with /z = 5, 6), as a function of time, while the animal 
moves with constant running speed. Black curves shows 
the thcta rhythm computed as the mean field of all the 
excitatory units of the network. 

The simulation is done at a constant animal speed v = 
1/ (400ms) (note that space is adimensional in our frame- 
work, and a unitary distance correspond to the distance 
between the centers of two adjacent receptive fields). The 
oscillatory input on each unit j is Ij{t) = Lj{t)cos{u)ot), 
where the amplitude Lj{t) is shown in fig. 2B as a func- 
tion of time. The figure shows that, for each place cell, 
the oscillation amplitude of the firing rate inc;reases as the 
animal approaches the center of its place field and de- 
crease as the animal left the center (firing tuning curve), 
while the place cell phase advanced at consistently earlier 
phases of the theta cycle as the animal pass through the 
cell's receptive field (phase precession). To quantify the 
phase precession phenomena, we compute the phase shift 
between the theta rhythm and the place cell activity in 
each cycle, and we plot in fig. 2 the theta phase as a func- 
tion of the animal position during running, together with 
the amplitude modulation. 

The model dynamics is in agreement with the exper- 
imental evidence of theta phase precession in two main 
points: 1) the total phase precession is less than 360" 
(O'Keefe Recce 1993, Skaggs 1996); 2) as the animal enters 
the receptive field the absolute phase of the cell activity 
with respect to the phase of maximal excitatory units pop- 
ulation activity (theta rhythm) is slightly less than 180°, 
in agreement with experimental finding of Skaggs et al 
(1996). 

Effects on nonlinearity in our general model have been 
analyzed in Scarpetta et al (2002), where we distinguished 
two class of nonlinearity. For the nonlinearity of class I (as 
the sigmoidal function shown in fig. 3B) we expect that 
the nonlinearity do not change critically the linear results. 
Indeed comparing the linear case in fig. 2 with the non- 
linear simulations results shown fig. 3 we see that the dy- 
namics is affected by the nonlinearity at large amplitudes 
but the phase precession is very well preserved. 

To check if in our model the theta phase is better cor- 
related with position than with time, we simulate the net- 
work dynamics with different values of the animal running 
speed. In fig. 4 the theta phase is shown as a function 
of animal position when animal runs with three different 
speeds {v = 1/800tos^^ black squares, v = l/400ms~^ 
red triangles, and v = l/200ms~^ green circles). As we 
expect, the theta phase is a unique function of the animal 
position, almost insensitive to the running speed, as long 
as the reaction time of the network is short compared with 
the time characteristic of the animal movement. While, 
when the speed is large the inertia-like effects becomes 
significant and a (small) dependence from the speed ap- 
pears (green circles follow a line with a higher slope) in 
the second half of the receptive field. 

From figure 4 we see that, when there's variability in 
running speed, the phase shows a better correlation with 



position than with time. This is in agreement with exper- 
imental results OS O'Keefe Recce (1993) and Huxter et al 
(2003). 

Finally, in order to test our model, we compare it with 
the recent experiments of Zugaro et al 2005. We have sim- 
ulated the effect of a perturbation that silences all units 
(both excitatory and inhibitory ones) and reset the phase 
of theta rhythm. The activity of two place cells, before 
during and after the network perturbation, is shown in 
fig. 5. The perturbation last for 200 ms. Despite the 
theta phase reset and the transient interruption of firing, 
the theta phase of the first cycle after the recovery is more 
advanced than the theta phase of the last cycle before the 
perturbation, as in the experiments (Zugaro et al 2005). 
The decrease of the peak firing rate is stronger in our sim- 
ulations than in the experimental findings of Zugaro et 
al (2005), that reported only a small but not-significant 
decrease in peak firing rate. To quantify the degree of 
conservation of the phase precession and the effect on the 
amplitudes of oscillation we shows in fig. 5 the plot of the 
theta phase as a function of the position, and the ampli- 
tude of the oscillation versus the animal position. Figure 
5C should be compared with fig. 2 in Zugaro et al (2005). 

Figs. 6 and 7 show the synaptic connections Jj/c and 
Wij learned in our model, as a function of pre-synaptic 
cell index k, for a post-synaptic place cells (i even) and 
an auxiliary units (i odd). As expected, the synaptic con- 
nections are translation invariant, so Jy depend only on 
the difference i — j and on the parity of i (even or odd) . 
We note that connections have both positive and negative 
values. Since they are supposed to be projections from ex- 
citatory pyramidal cells, the negative connections weights 
means effectively inhibitory connections, which could be 
implemented by additional inhibitory interneurons with 
very short membrane time constant. Note that the synap- 
tic weights are highly asymmetric. 

Finally, note that in our model the activity of all the 
units, both excitatory and inhibitory units, is modulated 
at the theta frequency. The network shows resonance in 
the theta range, and the inhibitory units, similarly to the 
place cells, exhibit phase precession (not shown), while the 
auxiliary units have a very weak place field with very little 
phase precession (not shown). 

5 Conclusions and Discussion 

In this paper we propose a theoretical model for theta 
phase precession in the hippocampus by applying a frame- 
work (Scarpetta Zhaoping Hertz 2002) that is a general- 
ization of the Hopfield model to oscillatory patterns. The 
learning rule is a generalization of the Hebb prescription, 
inspired to the STDP synaptic plasticity. 

The associative memory behaviour of our model, allows 
the network to show both amplitude and phase coded in- 
formations in response to inputs that carry information 
about current location only in the amplitudes modula- 
tion. 

Our model dynamics is in agreement with the experi- 
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mental evidence of theta phase precession in 3 main points: 

1) the model accounts for the observations that after 
a perturbation that transiently turned off neuronal dis- 
charges for a couple of theta cycles and reset the phase 
of the theta oscillations, the phase precession continued 
(Zugaro et al. 2005). This happens in our model, even 
if the external input to the network is silenced during the 
perturbation. This comes naturally from the fact that in 
our model, the oscillation phase of a place cell depends 
on both the network connections and the external input 
strength. Hence, when the external input strength rep- 
resents the spatial location (relative to the center of the 
place field), the oscillation phase of the place cell will, by 
the influence of the network connections, code the location 
relative to the center of the place field. 

2) In our model the theta phase depends on current po- 
sition of the animal, and is only slightly sensitive to the 
animal speed (almost insensitive at low speeds). There- 
fore in our model the theta phase is better correlated 
with position than with time spent in the place field, in 
agreement with experimental data (O'Keefe Recce 1993, 
Huxter et al 2003). This is achieved without need of any 
speed-dependent input, in contrast with the recent model 
of phase precession (Lengyel et al 2003) that assumes a dif- 
ferential input proportional to the instantaneous velocity 
of the rat. 

3) Both the initial theta phase as the animal enters the 
place field and the total phase precession generated by our 
network dynamics, are in quantitative agreement with ex- 
perimental evidence. In our model the total phase preces- 
sion is less than 360°, and the absolute phase of the cell 
activity w.r.t. the average pyramidal cells activity, when 
the animal enters the place field, is slightly less than 180°, 
in agreement with experiments of Skaggs et al. (1996). 
Note that our definition of the theta rhythm as the av- 
erage activity of all the excitatory units is in agreement 
with the definition used by Skaggs et al (1996). 

Finally note also that our model predicts both a mono- 
tonic phase precession and a unimodal firing tuning curve 
with a large amplification factor when the animal is in the 
center of the place field. This is in contrast with some 
other models. For example, models in which both fir- 
ing rate and phase are linear or monotonic functions of 
the same parameter (Booth and Bose 2001; Magee 2001) 
produce a monotonically increasing firing rate inside the 
place field if monotonic phase precession is preserved, or 
non-monotonic phase precession if unimodality of the fir- 
ing profile is to be preserved. 

Although there is some controversy as to whether theta 
phase precession is linear (O'Keefe Recce 1993) or if theta 
phase is an accelerating function of position (Skaggs et al 
1996, Yamagushi et al 2002), and if it's smooth or clus- 
terizcd, its trend is undoubtedly monotonic (O'Keefe and 
Recce 1993; Skaggs et al 1996, Yamagushi 2002), since 
there is a systematic progressive phase retardation as the 
animal passes through the place field. 

In contrast, experimental data indicates that firing 
rate changed in a non-monotonic waxing-waning manner 
(O'Keefe Recce 1993, Skaggs et al. 1996), even though 



there is some controversy as to whether phase precession 
accompanies skewed firing profiles (Mehta et al. 2002; 
Harris et al. 2002) or not (Huxter et al. 2003). 

Our model shares some similarities to the one of 
Tsodyks et al 1996, since both are based on the associative 
memory properties of the network. In contrast with the 
Tsodyks et al (1996) work, here the couplings are com- 
puted using a mathematical learning rule and the network 
dynamics is studied analytically in a framework (Scarpetta 
et al 2002) that is the generalization of the Hopfield model 
to the case of not-static but oscillatory patterns. 

A point of our model that still needs investigation is the 
origin of the imprinted patterns used during the learning 
mode. Of course, if there is no phase differences between 
units in the imprinted patterns, there would be no phase 
precession despite of learning and network dynamics. We 
are currently investigating how such phase differences arise 
in the imprinted patterns during learning, i.e. if they are 
generated by an initial asymmetry of the connections that 
produces the phase shift among units when the animal is 
moving, or if they arc inherited from earlier stages input 
during learning. Anyway, what we have proven here is 
that when the network is equipped with such a synaptic 
connectivity the network dynamics shows phase precession 
and phase-coded response to sensory input that do not 
have any phase difference between units activities. 

It has previously been suggested that phase coding oc- 
curs in the sensory input to the place cells (Burgess et al. 
1994), and it has been shown that under that model in- 
dividual place cells will show phase precession versus the 
average place cell activity (Burgess et al. 1993), as in our 
model. A difference with their model is that, even if we use 
a phase-coded activity during learning, once the network 
is equipped with the proper synaptic couplings, the phase 
precession occurs, when the animal run, without need of 
any phase coding in the sensory input. 
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Figure 1: A. Example of a learning kernel A{t) that satisfy 
the resonance conditions 18, at frequency lvq = lOHz. 
(Parameters 7 = /3 = 0.2rns~^, a = 0.14rns~^). The real 
(B) and the imaginary (C) part of its Fourier transform. 
In all the simulations we used the learning kernel shown 
in A. 



7 



time (ms) 







' i' 













If 





8000 

time (ms) 



B 



C 



Figure 2: Simulation of the network dynamics when ani- 
mal is moving continuously in time with constant running 
speed v=l/(200 ms). A: The two subplots show the activ- 
ities (in dotted red) of two place cell units {uj, j = 10, 12, 
whose place fields are centered at 5 and 6 respectively) as 
a function of time. Black curves show the ongoing theta 
rhythm activity (computed as the mean field of the ex- 
citatory activity of all the network). The external input 
to the units is Ij{t) = Lj(t)cos{ujQt), where the amplitude 
Lj(t), shown in B, is maximum when the animal is in the 
center of the receptive field of that unit j. Note that each 
place cell unit shows phase precession, i.e., relative phase 
shift between u{t) (red) and the theta activity (black), as 
the animal enters and then leaves the place field of this 
unit. Same behaviour for all other place cells not shown. 
C: Theta phase as a function of the current position of the 
animal, and amplitude of oscillation as a function of cur- 
rent animal position, for the place cell centered in location 
5. Theta phase is a bit less than 180° when the animal 
enters the receptive field, and the total phase precession 
is less than 360°. 
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Figure 3: Effects of nonlinearity in the phase precession 
phenomena dynamics. A. Simulation of the network dy- 
namics using the nonlinear activation function showed in 
B, when animal is moving continuously in time with con- 
stant velocity v=l/(200 ms) as in previous figure. Dotted 
red curve shows the place cell activity as a function of 
time. Black curves show the theta rhythm. Results are 
affected by nonlinearity mainly at large amplitude, but 
the theta phase precession is preserved, as shown in figure 
C showing the theta phase as a function of position, and 
amplitude of oscillation as a function of position for the 
place cell with receptive field centered in 5. 
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Figure 4: Theta phase of the place cell as a function of the 
animal position (A), and as a function of the time spent 
in the place field (B) at three different animal running 
speeds. During the place field traversal, for each cycle, we 
compute the phase of the place cell with respect to the 
theta rhythm. In each of the three simulations, the ani- 
mal is moving continuously in time with constant speed. 
The green circles correspond to speed v= 1/200 ms~^, the 
red tringles to v= 1/400 ms~^, and the black squares to 
v=l/800 ms~^ . Figure shows that the theta phase is much 
better correlated with position than with time. 
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Figure 6: Jik connection after learning N=10 locations, 
as a function of k. On the left, i = 10, and on the right, 
i = 11. 




Figure 7: Wik connection after learning N=10 locations, 
as a function of k. On the left, i = 10, and on the right, 
i = 11. 
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Figure 5: Phase precession persist after transient intra-hippocampal perturbation. Despite the theta phase reset and 
the transient interruption of firing, the phases were still correlated with the spatial position of the animal immediately 
after the recovery, in agreement with experiments (Zugaro et al 2005). A. Simulation of the network dynamics when 
animal is moving continuously in time with constant velocity v= 1/(200 ms) but a perturbation that silences all units 
(both excitatory and inhibitory ones) and reset the phase of theta rhythm is applied for 200 ms. Black curves show 
the theta rhythm. The dotted red line shows the activity of the place cells before during and after the network 
perturbation (place cell us, centered in location 4, in the upper subplot, and place cell miO, centered in 5, in the lower 
subplot). The theta phase and the firing amplitude are shown as a function of animal position when perturbation was 
applied, and in control conditions (i.e. without any perturbation), for the place cell centered in 4 (B) and the one 
centered in position 5 (C). 
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